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I. INTRODUCTION 

The long-range interaction between an atom and a per- 
fectly conducting wall is governed by the Lennard-Jones 
PJ static image potential 

V(z) = -^, (1) 

where z is the distance between the atom and the wall. 
The coefficient C3 in Eq. is the expectation value of 
the operator 

^ {xiXj + ViVj + 2ziZj) , 

in the atomic ground state. Here, = (xi, yi, Zi) is 
the coordinate of the «-th atomic electron with respect 
to the nucleus. For an atom with a spherically symmet- 
ric ground state, one can replace C3 by the equivalent 
expression 

C 3 = ^(0|i? 2 |0), 

where R = X)i r *- The Lennard-Jones interaction con- 
stant is important in models accounting for the finite con- 
ductivity of the wall material by Bardeen and Mavroy- 
annis [3|. Additionally, the wall- atom- wall interaction 
constant for small wall separation distances is propor- 
tional to C 3 H0. 

Precise values of C3 for lithium were obtained by Yan 
and Drake [(| from an elaborate configuration interaction 
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(CI) calculation and confirmed by an independent calcu- 
lation by Yan et al. The CI value of C3 for lithium 
is in close agreement with the value inferred from a vari- 
ational calculation by King 0. These accurate values 
of C3 for lithium are about 2% smaller than the value 
obtained from a Hartree-Fock (HF) calculation. 

An accurate semi-empirical value of C3 for sodium was 
also obtained by Kharchenko et al. p| from an analysis 
of the S-\ sum rule: 

^ = f<* 2 > = ?^ 

where /„ is the oscillator strength of the transition from 
the ground state to an excited state n. The quantities 
Eq and E n are energies of the ground state and excited 
state, respectively. This value differs from the HF value 
of C3 by about 10%. The more elaborate calculations 
by Per Jonsson et al. improve the agreement between 
theoretical and semi-empirical values for sodium some- 
what. 

Third-order many-body perturbation theory calcula- 
tions of C3 for all alkali-metal atoms and all-order singles- 
doubles calculations of C3 for Li, Na, and K were given 
by Derevianko et al. |j • The all-order calculations for Li 
and Na were in close agreement with with other precise 
values. More recently, Derevianko et al. ^(| deduced ac- 
curate theoretical and semi-empirical values of C3 for all 
alkali- metal atoms from oscillator-strength sum rules. 

In the present work, we use finite-field many-body 
methods to obtain values of (R 2 ) and make compar- 
isons with previous work. One advantage of the finite 
field-method is that all-order random-phase approxima- 
tion (RPA) corrections are included from the start. A 
second advantage is that n-th order corrections to the en- 
ergy give many-body perturbation theory (MBPT) cor- 
rections to matrix elements of R 2 normally associated 
with order (n+1). A third advantage is that matrix ele- 
ments of one- and two-particle operators are essentially 
trivial to obtain, in contrast to the lengthy calculations 



2 



ordinarily required. 



II. METHOD 

The method used here to evaluate the expectation 
value of the operator 



R 2 



N 



E r * = E r * +2 E r '-^ 

i i^j 



(2) 



is known as the "finite-field" method and widely used in 
quantum chemistry. We evaluate the ground-state expec- 
tation value of the operator by adding the scalar operator 
XR 2 to the many-electron Hamiltonian and calculating 
the resulting energy E(X). The desired expectation value 
is then given by 

/ , dE 

(R 2 ) = lim — . 

To evaluate the energy, we use standard many-body 
methods. 

The modified many-electron Hamiltonian may be writ- 
ten 



Mi 



i 



2A Tt • r, 



where, 



h a = ca-p + (3mc 2 + V nuc (r) + Xr 2 + U(r). 

Expressing the Hamiltonian in second-quantized form 
and normally ordering with respect to the closed core, 
we find 



H = E Q + U ■ aW- + o E v ^ kl : a i 



aja z a fc : 



ijkl 



where 



E(^F-C/) y - :at 



(3) 



In the above, equations, 

(Vhf). (J = E { V W« ~ u 



laaj J j 



with 



Vijkl 



'J 



— + 2Aiyr, 

7*12 



kl 



The lowest approximation to the wave function for a 
closed-shell atom is 

*0 = 4°6"-°n|0)- 



The expectation of the Hamiltonian in this state, which 
gives the first approximation to the ground-state energy, 
is 

<*„|ff |*0> = E = ?a + E f^HF - U 



In particular, if we choose U = Vhf, then the correspond- 
ing value of the closed-shell energy is 



E 



HF 



E 



i 



HF) 



III. RESULTS 



A. HF-level 



We consider an atom with one electron beyond closed 
shells. The core and valence energies are given at the HF 
level of approximation by 



£(HF) = £ 
a 

£l HF) = e,„ 



1 



(4) 
(5) 



where the summation index a ranges over closed shells, 
and where e v is the eigenvalue of the valence "frozen core" 
(Dirac)-Hartree-Fock equation. 

In setting up the HF equations, we add Xr 2 to the 
nuclear potential and add 2A r j • Yj to the two-particle in- 
teraction that defines the HF potential. The modified 
electron-electron interaction becomes 

1 °° r L 

— + 2An- r 2 = V 4TTPi(cos6') + 2Ar 1 r 2 Pi(cos6l). 

It follows that only the L — 1 term in the HF potential 
is modified and this term becomes 

vi(ab, r) — > vi(ab, r) 

/>oo 

+ 2Xr / r' [P a (r')P b (r') + Q a (r')Q b (r')] dr', 
Jo 

where P a {r) and Q a (r) are large and small component 
radial Dirac wave functions, respectively. As a practical 
matter, we choose A <C 1 for neutral atoms to main- 
tain some resemblance to the usual HF picture. The HF 
energy includes first-order MBPT corrections, together 
with all second- and higher-order RPA corrections. In 
columns 2 and 3 of Table [I] we list HF valence energies 
Ei HF ^ (A) and the HF core energies Ei HF ^ (A) as functions 
of A for the alkali-metal atoms from Li to Fr. 



B. 2nd-order MBPT 

We can easily go beyond the HF approximation and 
include the second-order MBPT corrections to the core 



3 



energy and to the valence energy 



e9 



^mnab^abmn 



ibrmi 



(6) 



E 



^mnvb^vbrnn 



Cb 



^mvab^abrnv 



abrn 



£6 



(7) 



In the above equations, indices a and b refer to core or- 
bitals, indices m and n refer to virtual orbitals, and v 
refers to the valence orbital. The second-order energies 
include those corrections to the matrix element of R 2 
usually associated with third-order MBPT - one interac- 
tion with R 2 and two Coulomb interactions. In columns 
4 and 5 of Table|I| we list E^ (A) and E {2) (A) for various 
values of A. 



C. Differentiation Formulas 

The energies are given on the grid 

A„ = (-3/i, -2h, -h, 0, h, 2h, 3h) 

with spacing h = 0.00002. To evaluate (-R 2 ), we make use 
of an hierarchy of successively more accurate Lagrangian 
differentiation formulas 



dE\ {3) 
dX ) 

dE\ {5) 
dX J 



dE\ 
~dX) 



(7) 



-2{E[2]-E[-2]) 
-108(^[2] - E[-2}) + 12(E[3] - E[-3]) 



to obtain dE/dX at A = 0. In the above, we designate 
E(X n ) by E[n]. The k-th approximation to the derivative 

(dE/dX)^ has an error proportional to h k . The first two 
of these formulas are given in Abramowitz and Stegun 
[H Chap. 25]. 

In Table [H] we show results of applying the differen- 
tiation formulas to the data in Table HJ The resulting 
values of (R 2 ) are numerically stable to about 4 digits 
for the cases considered. These values are compared with 
values from third-order MBPT and other accurate values 
in Table ITTT1 



D. Third-Order MBPT 

Expressions for third-order correlation corrections to 
core and valence energies of atoms with a single valence 



TABLE I: Variation of MBPT contributions to energies 
of alkali-metal atoms with the scaling parameter A, where 
H(X) = H + XR 2 . 



E. 



(HP) 



(A) 



E, 



(HP) 



(A) E { V 2) (X) Ei 2) (X) 



-0.00006 
-0.00004 
-0.00002 
0.00002 
0.00004 
0.00006 

-0.00006 
-0.00004 
-0.00002 
0.00002 
0.00004 
0.00006 

-0.00006 
-0.00004 
-0.00002 
0.00002 
0.00004 
0.00006 

-0.00006 
-0.00004 
-0.00002 
0.00002 
0.00004 
0.00006 

-0.00006 
-0.00004 
-0.00002 
0.00002 
0.00004 
0.00006 

-0.00006 
-0.00004 
-0.00002 
0.00002 
0.00004 
0.00006 



-0.1973918 
-0.1970331 
-0.1966760 
-0.1959663 
-0.1956137 
-0.1952625 

-0.1832811 
-0.1828627 
-0.1824466 
-0.1816209 
-0.1812111 
-0.1808033 

-0.1493919 
-0.1487522 
-0.1481185 
-0.1468683 
-0.1462511 
-0.1456389 

-0.1414238 
-0.1407050 
-0.1399943 
-0.1385951 
-0.1379058 
-0.1372228 

-0.1299267 
-0.1290615 
-0.1282090 
-0.1265379 
-0.1257176 
-0.1249066 

-0.1334806 
-0.1326676 
-0.1318658 
-0.1302922 
-0.1295191 
-0.1287542 



Li 

-7.2372588 
-7.2372410 
-7.2372232 
-7.2371876 
-7.2371698 
-7.2371520 
Na 

-161.8961446 
-161.8960647 
-161.8959849 
-161.8958252 
-161.8957454 
-161.8956655 
K 

-601.3789128 
-601.3786754 
-601.3784379 
-601.3779630 
-601.3777256 
-601.3774882 
Rb 

-2979.6664077 
-2979.6660442 
-2979.6656808 
-2979.6649539 
-2979.6645905 
-2979.6642272 
Cs 

-7786.6477893 
-7786.6472240 
-7786.6466589 
-7786.6455289 
-7786.6449640 
-7786.6443992 
Fr 

-24307.9714413 
-24307.9707289 
-24307.9700168 
-24307.9685926 
-24307.9678807 
-24307.9671690 



-0.0016257 
-0.0016327 
-0.0016397 
-0.0016535 
-0.0016603 
-0.0016671 

-0.0057754 
-0.0058073 
-0.0058389 
-0.0059014 
-0.0059323 
-0.0059630 

-0.0120214 
-0.0121561 
-0.0122889 
-0.0125494 
-0.0126774 
-0.0128041 

-0.0144304 
-0.0146265 
-0.0148194 
-0.0151967 
-0.0153817 
-0.0155647 

-0.0168211 
-0.0171307 
-0.0174334 
-0.0180225 
-0.0183103 
-0.0185943 

-0.0206103 
-0.0209570 
-0.0212968 
-0.0219595 
-0.0222839 
-0.0226043 



-0.0400988 
-0.0400995 
-0.0401001 
-0.0401014 
-0.0401020 
-0.0401026 

-0.3836680 
-0.3836734 
-0.3836788 
-0.3836896 
-0.3836950 
-0.3837004 

-0.7237274 
-0.7237770 
-0.7238267 
-0.7239262 
-0.7239760 
-0.7240257 

-1.8931869 
-1.8932719 
-1.8933569 
-1.8935271 
-1.8936122 
-1.8936974 

-3.1079015 
-3.1080755 
-3.1082495 
-3.1085978 
-3.1087722 
-3.1089466 

-5.8273192 
-5.8275533 
-5.8277875 
-5.8282563 
-5.8284909 
-5.8287257 



electron were given in |l2j | and applied to study ground- 
state removal energies of Cs and Tl in [l3j . In the present 
applications, these formulas are used to evaluate, effec- 
tively, fourth-order corrections to matrix elements of R 2 . 
Although we do not expect the third-order calculations 
presented in this section to be as accurate as the singles- 
doubles (SD) calculations given in the following section, 
it is in any case necessary to carry out third-order energy 

(3) 

calculations to determine E^ x t ra , the correction to the SD 
energies that accounts approximately for omitted triple 
excitations in the SD equations. 

Third-order corrections for lithium are dE y v °> /dX = 



4 



TABLE II: Values of {dE/dX)x=o = (R 2 ) for alkali-metal 
atoms as order of differentiation formula is increased. Step 
size: h = 0.00002. 



Order dEl HF) /d\ 



dE, 



(HP) , 



[dX dE { v 2) /dX dE { c 2) /d\ (R 2 ) 





T : 

Ll 








1 •? T A 1 O 


a o aa a 
0.8904 


A O A A C 

-0.3445 


A AO 1 E 

-0.0315 


18.256 


L(.( 410 


U.S9U4 


A 9 A A C 

-U.3445 


A A0 1 C 

-U.U315 


18.250 


1 T T A 1 K 

1 /. /41o 


A O AA A 
0.8904 


A O /I /I E 

-0.3445 


A A0 1 E 

-0.0315 


i o oct; 
18.256 




NT , 

IN a 








zU.04,30 


O AAOO 

3.9922 


-1.5631 


A oc; A A 

-0.2094 


OO O AO 

22.803 


on ^ /i ot 
z(J.04z f 


O AAOO 

3.9922 


-1.5629 


a nan a 
-0.2094 


OO O AO 

22.803 


on a a 07 
zU.04z f 


Q AAOO 

3.9922 


-1.5029 


A OfiA/1 

-U.2094 


OO QAO 

22.8U3 


31.2556 


K 

11.8718 


-6.5118 


-2.4864 


34.129 


oi o E O 1 

ol.zool 


1 1 C71 O 

11.8/ 18 


-0.5105 


O A Qd A 

-2.4864 


OA 1 OO 

34.128 


oi o COO 

ol.zooz 


1 1 C71 O 

11.8/ 18 


C K 1 A C 

-D.5105 


O A Qd A 

-2.4864 


OA 1 OO 

34.128 




Kb 








34.9792 


18.1725 


-9.4361 


-4.2549 


39.461 


34.9755 


18.1729 


-9.4336 


-4.2549 


39.460 


34.9755 


18.1730 


-9.4337 


-4.2549 


39.460 




Cs 








41.7773 


28.2500 


-14.7337 


-8.7108 


46.583 


41.7703 


28.2499 


-14.7276 


-8.7108 


46.582 


41.7704 


28.2498 


-14.7278 


-8.7108 


46.582 




Fr 








39.3401 


35.6056 


-16.5691 


-11.7207 


46.656 


39.3343 


35.6068 


-16.5634 


-11.7207 


46.657 


39.3344 


35.6073 


-16.5635 


-11.7207 


46.658 



TABLE III: Comparison of the present second-order finite- 
field (FF (2) ) values of R 2 with third-order MBPT values from 
Q and with semi-empirical (SE) values from [Tfl ]. The value 
for Li reported under SE is obtained by rounding the "exact" 
value 18.216. . . given in Refs. [6j,|7J. 



Element 



FF (2) 



MBPT 



SE 



Li 

Na 

K 

Rb 

Cs 



18.26 
22.80 
34.13 
39.46 
46.59 



18.26 
22.79 
34.05 
39.37 
46.35 



18.22 
22.45 
34.52 
40.92 
50.96 



-0.0297 and dE [ c ] /d\ = -0.0011. Adding these val- 
ues to the earlier second-order result, leads to (R 2 ) = 
18.2250 for lithium. This slightly improves the agreement 
of MBPT with the exact nonrelativistic value. However, 
better agreement can be achieved in the SD approxi- 
mation. Therefore, here we only calculate dE) ;y : tla / dX. 
These contributions are 0.00252 for Li, 0.00201 for Na, 
0.72532 for K, 1.10114 for Rb, 2.32450 for Cs, and 
2.59698 for Fr. 



E. All-Order Singles-Doubles 

The singles-doubles (SD) equations, also referred to 
as all-order pair equations [141 Chap. 15], provide a 



method of including important correlation corrections to 
the atomic wave function to all-orders in perturbation 
theory. One solves a set of coupled equations for single 
excitation coefficients p ma , p m v and double excitation co- 
efficients pmnab, Pmnva of the HF ground state [see HH 
for example] . Once these expansion coefficients have been 
determined, the correlation correction to the core energy 
AE C is given by 



AE r 



fcmn Pmnab- 



(8) 



and the correlation correction to the valence energy is 
given by 



i Pma 



ma 

H - ^ ^ ^abvmPmvab H~ ^ ^ ^vbmnPri 



(9) 



The core energy is exact through third order in MBPT 
and contains important fourth- and higher-order correc- 
tions. The valence energy also includes important fourth- 
and higher-order corrections but is missing small third- 
order corrections [written out explicitly in IpoT ] referred 



to as E, 



(3) 



These missing terms have their origin in 



omitted triple excitations (single-valence - double-core 
excitations) of the HF ground state. Numerical values of 

(3) 

-^cxtra f° r the alkali-metal atoms are given at the end of 
Section ITTTnl 

a. Lithium: Calculations for Li include all partial 
waves with / < 6. To estimate higher I contributions we 
use Aitken's 5 2 method. Table IIVI shows contributions 
to (R 2 ) evaluated with l max ranging from 2 to 5. The 
final extrapolated value (i? 2 ) = 18.213 from Table |IV| 
differs from the "exact" nonrelativistic value (18.216004) 
for lithium given by Yan and Drake but is in precise 
agreement with an earlier SD result by Derevianko et al. 
|9j. The small difference with the exact nonrelativistic 

(3) 

value is dominated by the contribution from E^ x [ la , eval- 
uated in the previous subsection, which has the value 
dE^^/dX = 0.0025. When this correction is added to 
the SD result 18.2130 for lithium, we obtain the value 
18.2155, differing from the exact nonrelativistic result by 
only -0.0005. The residual difference has the sign and 
order of magnitude expected for a relativistic correction 
to R 2 . 

b. Other alkalis In TableEl we show the derivatives 
of valence and core energies of alkali atoms from Li to 
Fr calculated in the SD approximation with Z max = 6 
as the order of the differentiation is increased. We also 
include the contribution from the missing third-order en- 

(3) 

ergy £^xtra evaluated in the previous section. 

The SD result for sodium (R 2 ) = 22.6425(3) agrees 
well with the earlier SD result 22.6293 from Derevianko 
et al. ;9| and with the semi-empirical value 22.65 from 
Kharchenko et al. 0. Note however, that present results 
for all alkali atoms other than lithium are substantially 



5 



TABLE IV: Calculated values of (R 2 ) for Li as the number of 
partial waves £ m ax included in the SD equations is increased 
are tabulated along with extrapolated values obtained by ap- 
plying Aitken's 5 2 method to i max = (2, 3, 4) and (3, 4, 5). 



^max 


dAE v /d\ 


dAEc/d\ 


(R 2 ) 


2 


-0.37905 


-0.03284 


18.2200 


3 


-0.38362 


-0.03279 


18.2155 


4 


-0.38509 


-0.03279 


18.2140 


5 


-0.38570 


-0.03279 


18.2134 


2-3-4 


-0.38579 


-0.03279 


18.2133 


3-4-5 


-0.38612 


-0.03279 


18.2130 



larger than semi-empirical values obtained in Ref. [To| 
(see Table HTU. 

The resulting values (R 2 ) from the SD calculation, 
which are our most accurate predictions, are summarized 
in Table EU 



Conclusion 



TABLE V: Values of (dE/dX) x=0 in the SD+E^ ra approx- 
imation as the order of differentiation formula is increased. 
Step size: h — 0.00002. The SD equations included all partial 
waves with / < 6 for Li, Na, K, Rb, Cs and with / < 5 for Fr. 



Order 


dE< HF > 


dE< HF > 


dAE v 


dAE c 




(R 2 ) 


dX 


dX 


dX 


dX 


d.X 








Li 








3 


17.7418 


0.8904 


-0.3860 


-0.0328 


0.0025 


18.216 


5 


17.7415 


0.8904 


-0.3859 


-0.0328 


0.0025 


18.216 


7 


17.7415 


0.8904 


-0.3859 
Na 


-0.0328 


0.0025 


18.216 


3 


20.6433 


3.9922 


-1.7038 


-0.2922 


0.0020 


22.642 


5 


20.6427 


3.9922 


-1.7037 


-0.2922 


0.0020 


22.641 


7 


20.6427 


3.9922 


-1.7037 
K 

-5.9711 


-0.2922 


0.0020 


22.641 


3 


31.2556 


11.8719 


-2.2866 


0.7255 


35.595 


5 


31.2531 


11.8718 


-5.9703 


-2.2866 


0.7253 


35.593 


7 


31.2532 


11.8718 


-5.9703 
Rb 


-2.2866 


0.7253 


35.593 


3 


34.9792 


18.1725 


-8.2136 


-3.6888 


1.1014 


42.351 


5 


34.9755 


18.1729 


-8.2119 


-3.6888 


1.1011 


42.349 


7 


34.9755 


18.1730 


-8.2117 

Cs 


-3.6888 


1.1011 


42.349 


3 


41.7773 


28.2500 


-11.6972 


-6.6590 


2.3254 


53.997 


5 


41.7703 


28.2499 


-11.6942 


-6.6590 


2.3245 


53.991 


7 


41.7704 


28.2498 


-11.6943 
Fr 


-6.6590 


2.3245 


53.991 


3 


39.3401 


35.6056 


-12.4105 


-8.5968 


2.5970 


56.536 



In this paper we present the most complete fully ab ini- 
tio all-order calculations of the Lennard-Jones interaction 
constant C3 for alkali-metal atoms. Incorporating of the 
rescaled R 2 operator into original Hartree-Fock Hamilto- 
nian allows us to stay within standard SD technique while 



TABLE VI: Final results for (R 2 ) for alkali atoms. 



Li 


Na 


K 


Rb 


Cs 


Fr 


18.216 


22.641 


35.593 


42.349 


53.991 


56.536 



also including important subclasses of higher-order con- 
tributions. Results for Li agreeprecisely with the "exact" 
CI results of Yan and Drake [g , while results for other 
alkali atoms are probably the most accurate available to 
date. 
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